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I. INTRODUCTION 



To describe phenomena like heavy ion collisions or the electroweak phase transition in 
the Early Universe, the methods of quantum field theory must be applied to finite densities 
and temperatures. The most reliable approach, especially for the case of QCD, are lattice 
calculations in Euclidean space-time (for a recent overview see e.g. JIJ). Unfortunately the 
results are only applicable to equilibrium processes. Dynamical processes in quantum field 
theory, on the other hand, are most commonly studied using perturbation theory. In the- 
ories with interacting massless degrees of freedom, however, naive perturbation theory in 
powers of the coupling constant breaks down at high temperature due to infrared diver- 
gences caused by the exchange of massless modes. To obtain finite results one must resum 
classes of Feynman diagrams which lead to dynamical mass generation and thus regulate 
the infrared singularities. To do this systematically one must disentangle the influence of 
different momentum scales; in QCD at finite temperature, these are characterized by T, gT, 
and g 2 T where the last two give the order of magnitude of the dynamical masses for the 
originally massless electric and magnetic gluonic modes, respectively. An improved one-loop 
perturbation scheme based on these principles, called hard thermal loop resummation, was 
developed by Braaten and Pisarski || (see also |J for an overview). Using effective field 
theory methods, the above separation of the different momentum scales was achieved by 
Braaten and Nieto for static quantities at high temperature Q. 

At high temperature, it is reasonable to expect that phenomena which are dominated 
by long wavelength modes can be reliably described by classical field theory. This is the 
basic motivation for the work presented here where we study the classical dynamics of the 
coupled Yang-Mills-Higgs system. It is well known that the classical non-Abelian gauge field 
equations are nonintegrable and exhibit dynamical chaos ||. Therefore one has to rely on 
numerical methods to study the evolution of classical fields. 

The real time evolution of classical SU(iV) Yang-Mills fields on a lattice, for N = 2,3, 
was recently studied in Ref. ||. Starting with randomly chosen initial conditions for the 



2 



gauge fields on each lattice link it was found that the system shows chaotic behavior, i.e. 
the distance in phase space between two adjacent configurations increases exponentially in 
time. For systems with sufficiently weak coupling, the time-averaged growth rate (maximal 
Lyapunov exponent) An was shown to scale linearly with the total energy of the system. 
The slope depends on the number N. Even if one starts with highly off-equilibrium config- 
urations, the gauge fields thermalize very fast, i.e. the energy distribution over the lattice 
soon exhibits a thermal shape. Therefore the total energy of the system can be related to 
the associated temperature, and the scaling of the maximal Lyapunov exponent with the 
total energy can be reinterpreted as a scaling with the temperature T - at least in the weak 
coupling regime. Precisely the relation 

Ano ~ g 2 Ta (1.1) 

which connects dimensionless quantities was established. Here a is the lattice spacing and 
g the coupling constant. One surmises that the actual value of the proportionality factor 

An« n n 

c N = -— - (1.2) 
9 T a 

survives in the continuum limit a — > 0. 

For the classical SU(2) and SU(3) gauge field systems, it was found in || that the 
maximal Lyapunov exponent coincides numerically with twice the damping rate 7n of a 
gauge boson at rest as calculated quantum field theoretically 0] using thermal perturbation 
theory with hard thermal loop resummation: 

An ° 2^. (1.3) 



g 2 Ta g 2 T 

At first sight it is surprising that the Lyapunov exponent which is deduced from the classical 
field dynamics by an exact (i.e. non-perturbative) numerical simulation should be related 
at all with the plasmon damping rate which is calculated perturbatively in quantum field 
theory. However, the close relation between linear response functions calculated in the hard 
thermal loop approximation and classical transport phenomena has been known for many 



years |7]||, indicating that the hard thermal loops are essentially classical (see also 0). In 
addition it was argued in |TIJ that the damping of a gluon and the exponential growth of 
the distance between two adjacent classical trajectories are related. So far there is, however, 
no rigorous proof for the relation ( |1.3| ). 

An empirical approach to the question whether the coincidence between the Lyapunov 
exponent and (twice) the damping rate is only accidental would be to simulate other systems 
on the lattice and also calculate perturbatively their damping rates. A natural extension 
of the pure SU(iV) system is to couple it to a matter field. In this article we study the 
real-time evolution of the coupled SU(2) Yang-Mills-Higgs system. Again we are interested 
in the thermalization properties and the chaoticity of the fields. To compare with the high 
temperature results from hard thermal loop calculations we concentrate on the symmetric 
phase, i.e. we simulate the Higgs behavior in the electroweak plasma above the phase tran- 
sition in the 8 W = limit. Our work was partly motivated by the recent study presented in 



Ref. [DJ in which the question of the connection between the damping rate and the Lya- 



punov exponent in such a system was also addressed [II]. This paper provides in particular 
the analytical results for the plasmon damping rates for the gluon and Higgs fields in the 



hard thermal loop approximation. However, in Ref. [IT] the thermalization properties of 
the classical fields were not taken into account, and thus the determination of the temper- 
ature of the system was not dealt with in a satisfactory manner. The correct extraction of 
the temperature is, of course, crucial to establish relations like ( |1.1|) . We postpone a more 



detailed assessment of the conclusions of Ref. |TT[ until after presenting our own results. 

The organization of the paper is as follows: In Sec. |H| we formulate the problem. The 
thermalization properties of the fields are studied in Sec. |T|. In Sec. [IV] we compute the 
Lyapunov exponents, and we discuss our results and compare them with previous work in 
Sec. [VJ 
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II. FORMULATION OF THE PROBLEM 



In this Section we present the realization of the coupled SU(2) Yang-Mills-Higgs system 
on the lattice, describe the initialization of the field variables, discuss the allowed ranges of 
the relevant parameters in the weak coupling and continuum limits, define the Lyapunov 
exponents and describe how we measure them. 



A. The Yang-Mills-Higgs system on the lattice 



In the temporal gauge Aq = 0, the Hamiltonian of the coupled SU(2) Yang-Mills-Higgs 
system in the continuum is given by 



H — Jd 3 x [lEfE? + \B?B$ + + (A$) + (A$) + A($ f $) 2 ] 



(2.1^ 



The dot denotes the time derivative, and 



V 



\ 

) 



(2.2) 



>)2 - i<j>3 

is a charged Higgs field in the fundamental representation of SU(2). The gauge field is 
described by vector potentials A$(x). Ef = —Af denote the electric fields, Bf the magnetic 
fields, and A the covariant derivatives. 

As mentioned before, the equations of motion derived from the Hamiltonian ( |2.1| ) are not 
integrable. One has to rely on numerical methods to study the time evolution of the fields. 
For numerical purposes we discretize the spatial coordinates by using a cubic lattice with 
N 3 sites. Time is kept as a continuous variable. For details about the numerical simulation 
of the real-time evolution of the system we refer to ||. 

To realize the coupled SU(2) Yang-Mills-Higgs system on the lattice without losing its 



gauge symmetry properties we use link variables [12,13] 



U Xji = exp (~igaA?(x) r ( 



(2.3) 
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instead of vector potentials A\. Here r c are the Pauli matrices. The covariant derivative 
acting on the Higgs field is replaced by the difference 

-(U x ,i$ x+i -$ x ) . (2.4) 

(X 

The electric field strength on the lattice is related to the link variables and their time 
derivatives by 

EU = -^{T a U x ,iUl,d- (2-5) 
ya 

The magnetic energy is expressed in terms of the so-called plaquette operator U X;i j, which 
is the product of all four link variables on an elementary plaquette with sites (x, x + i, x + 
i + j,x + j): 

with U x -i = Ul_ii. The links are directed and hence the plaquettes are oriented. In the 
continuum limit the plaquette variable U X jj is related to the local magnetic field B£ k by 

U x>ij = exp (-\iga\ ijk Bl k T«) . (2.7) 

It is useful to write the SU(2) matrices in the form 



U = u - ir a u a = 



( . \ 

u - tu 3 , —u 2 - lUi 



(2.8) 



u 2 - lUi , U + IU 3 J 

where the Ui are four real numbers, which can be thought of as components of a quaternion. 
Since det U — 1, each link variable U can be expressed in terms of three angles, pc, $q, ipc, 
characterizing a three dimensional hypersphere: 

U = cos — - ip a G r a sin ^ (2.9) 

with a three dimensional unit vector 

Pg = (sin-# G cos(^G ; sin-# G sin^ G , cos-^g) • (2-10) 
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Similarly the four independent components of the Higgs field can be represented as a quater- 
nion: 

\ 



$ = 0o - ir a (p a 



V 



(2.11) 



- #3 , ~</>2 - i(/>l 
02 - i(f>l , 00 + ^03 

In contrast to the link variable U the Higgs quaternion <3> has arbitrary length R defined by 



/ 



J R 2 = 02 + 02 + 2 + 2 = | tr ( $t $) 



(2.12) 



Again we can introduce three angles to get the representation 



$ = R (cos ^ - ip£r a sin ^ 



(2.13) 



Using the quaternion representations for the link variables and the Higgs fields, the Hamil- 
tonian of the coupled SU(2) Yang-Mills-Higgs system on the lattice reads 



H — H e \ + H mag + iJ kin + H pot , 



(2.14) 



where we define the electric energy 



(2.15) 



the magnetic energy 



(2.16) 



the Higgs kinetic energy 



Hkin = X! e kin(^) , 



(2.17) 



and the Higgs potential energy 



Hpot = e PotA x ) + e P ot,4(z) + e int(x, i) , 

x x x,i 

and the e's are energies associated with the respective site, link, or plaquette: 



(2.18) 



e el (*)=a 3 ]T±£ x y^, 

i 

g d 

e kin (x) =a 3 [|tr ($£$ x ) 
e P ot, 2 (x) = 6a [|tr($£$ x ) 
e P ot,4(x) = a 3 A [|tr 
eint(x,i) = -atr \^ x U x ^ xM 



(2.19a) 

(2.19b) 
(2.19c) 
(2.19d) 
(2.19e) 
(2.19f) 



To derive the last equation we have exploited the unitarity of the link variable U and the 
fact that the trace of a quaternion is real. 

The Hamiltonian Q2.14 ) is scale invariant. To see this we define the dimensionless vari- 
ables 



t = t/a, H = g 2 aH } E a xi = ga 2 E a xi , $ x = ag$ x , and X = X/g . 



In terms of these we get 



(2.20) 



'2.2V, 



+ E | tr ($l$ x ) + E tr - <I'ir,,,<I>,./ + A E I [tr 

where the dot now denotes the derivative with respect to the scaled time t. Therefore the 
only parameters in the system are the scaled energies and the scaled Higgs self coupling A 
while the lattice spacing and the gauge coupling constant do not appear explicitly. 

It is easy to derive the equations of motion from the lattice Hamiltonian ( |2.21| ) and the 
relations (12.51) and 



(2.22a) 
(2.22b) 
(2.22c) 

(2.22d) 



(u a U = -\ [E a x .(u ) Xji + e ahc E b x .(u c ) x , 



k = -A tr ('I-,«I» r ) 'I', • (i <!•, • E { U *,i ®x+i + / ; 



Here 

^r(<i.;.7"r,,,<i.,.,) (2.23) 

is the lattice version of the gauge-covariant isospin current of the Higgs field. 

In addition to satisfying these equations of motion the physical solutions must obey 
Gauss' lawQ 

DfE\ - \%g ($V a $ - $V$) = (2.24) 

with the covariant derivative Df b in the adjoint representation. The expression on the l.h.s. 
of (|2.24 ) is a conserved quantity since it commutes with the Hamiltonian of the system. 



Thus to satisfy Gauss' law for all times we must only make sure that our initial conditions 
obey ( pg) . 

For a study of the thermalization properties of the different degrees of freedom the energy 
transfer between different sectors of the system is an important quantity. As we will see 
below the energy transfer from the gauge part of the system to the Higgs part is of particular 
interest. It is given by 

j (H el + H mag ) = -j (H kin + H pot ) = £ /-,./;„ (2.25) 

x, i 

as can be easily derived from the equations of motion. The local Higgs isospin current j° Li 
is defined in (ggg ). 

B. Initialization, weak coupling and continuum limit 

We initialize the system in the same way as in by choosing random initial configura- 
tions: 



For brevity we present only the continuum version of Gauss' law. 
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The isospin directions of the Higgs fields, characterized by the angle variables pn, ^h, 
<Ph are chosen randomly for each lattice site. To keep control over the initial potential 
energy of the Higgs fields the length of each Higgs quaternion is not randomly chosen. 
Instead we take for simplicity the same length R for all & x initially. Obviously the 
initial Higgs potential energy ( [2.18 ) increases with R. 



The isospin directions of the gauge fields are also randomly chosen. In this case, 
however, this does not mean that all the angles are arbitrary. The reason is that 
not only the isospin direction of the vector potential A but also its amplitude enters 
the link variable U as an angle variable. It is easy to see that pc is connected to the 
amplitude while $g and <Pg characterize the isospin direction. Therefore the latter two 
angles are initialized with random values on each link. To limit the initial magnetic 
energy we restrict the gauge field amplitude. Thus we introduce a parameter 5 which 
controls the initial magnetic energy and choose pc o n each link randomly but within 
the range [0, 2irS]. For 5 = the initial magnetic energy also vanishes. For small S the 
relation i/ mag ~ 5 2 holds. 

To satisfy Gauss' law ( |2.24| ) in an easy way we choose initially vanishing field momenta 
on each link, i.e. E = = This implies that the initial electric energy and the Higgs 
kinetic energy, fl2.15|) and Q2.171 ) respectively, vanish. Note that with non- vanishing po- 



tential but vanishing kinetic energies the system is initially in a highly off-equilibrium 
state. 

When initialized in this way the system's further evolution is determined by the three 
parameters 5, R, and A, or by the initial magnetic energy and the initial Higgs potential 
energy (quadratic and quartic part), respectively. 

To make contact with perturbative calculations it is necessary to stay in the weak cou- 
pling regime. Furthermore, weak coupling is required for the classical calculation to be a 
valid approximation to the full quantum treatment. This puts constraints on the range of 
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the physical parameters 5,R,\. The scale invariance of the Hamiltonian allows to fix the 
gauge coupling at g — 1 and the lattice spacing at a = 1. 

Weak coupling within the gauge sector is realized if the link variables are in the vicinity 
of unity (cf. ( |2.3| )). In the initial state this is achieved by choosing 5 small which implies that 
Pq in (|2.9|) is small. On the other hand, as first noted in ftM and recently worked out in W5 



one encounters finite size artifacts in the pure SU(2) system when using too small values of 
5 on a finite lattice. To explore the weak gauge coupling regime one is thus restricted to a 
window of not too small and not too large values for S, or one must take very large lattices. 
With the latter option one soon runs into computer limitations. 

From the Hamiltonian ( [2.21|) one sees that weak coupling within the Higgs sector is 



ensured if the quadratic terms dominate the quartic one. A system with weak Higgs self 
coupling can thus be initialized by requiring 

\R* < R 2 . (2.26) 

To ensure weak coupling between the gauge and the Higgs fields, the Higgs isospin current 
j in ( [2.22c ) should be small compared to the magnetic energy contribution^]. Estimating 



both contributions by the (maximal) field amplitudes we get the condition 

R 2 <5. (2.27) 

The continuum limit also puts constraints on the possible parameter ranges. First, as 
seen from Eq. ( |2.3| ), to recover the continuum limit from the lattice formulation one needs 
gaA(x) — » 0. If the lattice simulations are supposed to yield a reasonable approximation to 
continuum behaviour, 5 must therefore not be large. Another lattice artifact shows up if the 
total energy in the system is too high. Due to the compactness of the space of link variables 



the magnetic energy on the lattice ( |2.16| ) is limited from above by 



(^mag)max = 4^[l-|tr(-l)] =24iV 3 . (2.28) 



Tn the Abelian and continuum limits the first contribution on the r.h.s. of Eq. ( 2.22c ) is V x B. 
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If the total energy in the system is much larger than this value, equipartition of the energy 
among the magnetic, electric and Higgs degrees of freedom may thus not be possible. 

All these considerations must be taken into account when selecting the initial parameters 
S and R as well as the value A for the Higgs self-coupling. In our calculations we have made 
sure to choose these parameters in such a way that we are simultaneously in the weak 
coupling regime and close to the continuum limit. 



C. Lyapunov exponents 

In the next Section we will study the thermalization properties of systems initialized as 
described above. The self-thermalization of the fields is strongly influenced by the chaoticity 
in its dynamics which is characterized by a set of positive Lyapunov exponents. We will here 
concentrate on the largest Lyapunov exponent associated with the most unstable dynamical 
mode. In principle one can define two such exponents, Ac and Ah, for the gauge and 
Higgs fields, respectively, by using two different distance measures in the space of field 
configurations. The first measures the growth rate of the Euclidean distance in the gauge 
field sector, 

Dl[U, f/'] = EE [{{u*U - KM* + i«U - «Uf] , (2-29) 

x,i a=0 

where U and U' denote the gauge fields associated with two solutions of the equations of 
motion ( |2.22| ) resulting from two initially very close field configurations. The other Lyapunov 
exponent measures similarly the growth rate in the Higgs sector, 

DV& , = E E [((&«). - (€)*) 2 + ((*«). - (<^) 2 ] , (2-30) 

where $ and $' are the corresponding Higgs fields from the same two solutions. As we will 
see in the following, the two Lyapunov exponents Ac and Ah turn out to be identical. 
As usual the Lyapunov exponent is defined as 

A = lim lim -In (2.31) 

i-»D(0Ht D(0) K ' 
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To avoid the problem of distance saturation || in directly measuring the Lyapunov exponent, 
we use the rescaling method for the evaluation of the large-time limit required in ( [2.31| ) (see 
H for more details). 

In principle other definitions for the distance between field configurations are possible. 
In particular, one might wish to replace the Higgs distance ( |2.3U| ) by a gauge invariant 
definition like e.g. |6|.[TT 



XK-<1- (2.32) 

X 

However, the Euclidean distances defined above are more appropriate for the rescaling 
method. For the pure SU(2) gauge system it was carefully checked that different distance 
measures give rise to the same Lyapunov exponent |TH . 



III. THERM ALIZATION PROCESSES 

In this Section we investigate various thermalization processes in the coupled SU(2) Yang- 
Mills-Higgs system. We will encounter different types of thermalization processes which 
operate on different time scales. In particular we find that, starting with the initialization 
described in Sec. [TI], the gauge and the Higgs subsectors each thermalize rather rapidly, while 
the equipartition of energy between them requires a much longer time. In the following two 
subsections we first present results on the fast thermalization processes associated with the 
gauge and the Higgs subsystems, and then discuss the system's long-time behaviour. For 
reasons given in Sec. ||our study is confined to certain parameter ranges which are consistent 
with weak coupling and proximity to the continuum limit. 



A. Short-time thermalization and relaxation times 

Thermalization processes in a large system can be studied by following the evolution of 
the energy distribution over the microscopic degrees of freedom. In general, the distribution 
should approach that of a Gibbs ensemble upon equilibrium. As shown in the cases of pure 
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SU(2) and SU(3) gauge fields 0, by measuring appropriate energy distributions, one can 
determine a parameter T a which we call the apparent temperature henceforth. The apparent 
temperature T a is related to the system's true temperature T, but generally differs from T, 
if the degrees of freedom of the system are not all independent, such as for a gauge theory. 

We will test the degree of thermalization in various subsectors of the system by comparing 
the distributions for the associated energies over the lattice with a Boltzmann distribution 

P(e) = JV /(e)exp(-e/T a ). (3.1) 

Here e is the energy associated with the selected degrees of freedom (e.g. the magnetic or 

electric energy of the gauge fields), and /(e) is an appropriate phase space factor, specified 

below. No is a normalization constant. The distribution P(e) is obtained in the simulation 

by sampling, at a fixed time t, the variable e over all sites or plaquettes on the lattice; it 

gives the probability, averaged over the whole lattice, to find the value e at a given, though 

arbitrary, lattice point. If this distribution has the Boltzmann shape ( |3.1| ) the system has 

lost all initial order in the distribution of e, and the entropy has reached a maximum under 

the constraint of fixed total energy associated with the selected degrees of freedom. 

We will present measurements for four different energy distributions, the distribution 

P e i(e) for the electric energy e e i of the gauge fields, the distribution P mag (e) of the magnetic 

energy e mag , as well as the distributions Pkm(e) and P po t,2(e) of the kinetic and potential 

energies of the Higgs field, ekin and e pot> 2, respectively, as defined in Eqs. (|2.19|) . In the weak 

coupling limit the respective phase space factors /(e) are given by 

fe 2 forP cl (e), 




e forP kin (e), 
e for P po t i2 (e). 



A detailed calculation of the phase space factors can be found in Appendix [A]. 

By plotting ln[P(e)//(e)] against e one can read off the apparent temperature T a . It 
was discussed in PJT7[| that T a is not in all cases identical with the true temperature of the 
system. Writing 
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T a = C(T)T, 



(3-3) 



one finds 



f 1 for Pei(e) 



C 



C m for P mag (e), 
1 forP kin (e), 
k 1 for P P ot,2(e)- 



(3.4) 



C m is a number varying between 2/3 and 1 as T rises from low to high values ||. In Ref. it 
was argued that this can be understood in terms of the changing number of effective degrees 
of freedom which at low temperatures is only 2/3 of that at high temperatures, due to the 
increasing importance of longitudinal plasma modes at high temperature. More details on 
the relation between the apparent temperature and the actual temperature can be found in 
the Appendix [A]. 

By numerically evolving the initial configurations forward in time, we can thus follow 
the system's thermalization. At each moment, we can measure its energy distributions and 
then read off its characteristic temperatures from the measured distributions, using ( |3.3|) 
for correct normalization. In Fig. [1] several characteristic temperatures of the system are 
plotted as functions of time, for short times below t — 6 (in lattice units). To ensure weak 
coupling and closeness to the continuum limit we choose 5 = 0.2, R = 0.2, and A = 1. In 
this case it was found that C m ~ 0.687. 

The curves in Fig. [I] demonstrate a fast thermal equilibration within the gauge field 
sector between the electric (solid line) and magnetic (dashed line) energy distributions. 
After some transient oscillations they approach, on a typical time scale of 1 lattice unit, 
a common temperature T e \ w T mag w 0.47. On roughly the same time scale the kinetic 
(dotted line) and quadratic potential (dot-dashed line) parts of the Higgs field energies also 
equilibrate, but at a much lower temperature T kin ps T pot 2 ~ 0.07. This indicates that in 
the coupled SU(2) Higgs system the gauge and the Higgs subsector each undergo a separate, 
fast equilibration, starting from a highly off-equilibrium initial state. 
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The approximate equality of the two associated relaxation times, tq for the gauge sector 
and r H for the Higgs sector, can be easily understood if one makes the following two hy- 
potheses: 

1. Thermalization in a chaotic system is driven by the most unstable modes and thus dom- 
inated by its largest Lyapunov exponent. 

2. The gauge and the Higgs sectors are associated with the same largest Lyapunov expo- 
nent. 

While it is difficult to find a rigorous proof for the first hypothesis, such a behaviour is 
intuitively expected. The second hypothesis can be checked in our simulation and is found 



to be correct (see Sec. p"VP . 

At short times, each subsector settles down at a common temperature for all its degrees 
of freedom, but the temperatures of the two subsectors differ widely. This can be traced to 
the different initial energy content in the two subsectors. On average, the energy per degree 
of freedom in the gauge sector is initially 0.478, while that in the Higgs sector is only 0.060. 
These specific values are, of course, a result of the parameters chosen in our initialization 
procedure. They would be different if other parameter values were selected. 

The larger total gauge energy per gauge degree of freedom compared to the total Higgs 
energy per Higgs degree of freedom results in the higher temperature of the gauge sec- 
tor. Equilibration between the two subsectors apparently takes much longer. Since both 
subsectors are apparently equally "chaotic", the relative stability of the two temperatures 
of the subsystems indicates that energy exchange between the gauge and Higgs sectors is 
suppressed. This will be discussed quantitatively in the following subsection. 



B. Long-time equilibration and relaxation time 

As seen in Fig. [I], the coupled SU(2) Yang-Mills-Higgs system shows a very rapid ther- 
malization within each subsector, with the relaxation times being of the order of 1 lattice 
unit. However, it turns out that the mutual equilibration between the two sectors is rather 
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slow and requires much longer time. To study the long time behavior of the coupled system, 
we follow the evolution of the energies. 

In Fig. 0, we plot as a function of time (on a logarithmic scale) the following energies: 
E e \ (solid), E mag (dotted), Eq (short dash), E^l (long dash), and E t (dot-dash). Here E e \ is 
the electric energy per gauge degree of freedom: 

-Emag the magnetic energy per gauge degree of freedom: 

-Emag — g^y 3 H ma _ g , (3-6) 

En the total Higgs energy per Higgs degree of freedom: 



Eq the total gauge energy per gauge degree of freedom: 

E G = E e \ + E mSjg , (3.8) 

and Et the total energy per degree of freedom in the entire system: 

E t = (QE G + AE u )/10. (3.9) 

N 3 is the number of lattice sites. Note that each lattice site is associated with 10 degrees 
of freedom: 6 gauge degrees of freedom and 4 Higgs degrees of freedom. The parameters 
used in the calculation are the same as in Fig. [|. In the weak coupling regime, E^ is the 
temperature of the Higgs sector as Eq is that of the gauge sector, while E t will be the final 
temperature of the entire system after the coupled gauge and Higgs sectors reach complete 
equilibrium. 

Again, we notice that the energy transfer from the magnetic sector to the initially un- 
populated electric sector occurs on a time scale of the order of unity, which agrees with 
Fig. |l| The same has been seen for the Higgs kinetic energy and the Higgs potential energy 
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(not shown here). However, it takes a very long time for the initially hotter gauge sector to 
transfer some of its energy to the cooler Higgs sector. As shown in Fig. 0, Eq (short dash) 
and .Eh (long dash) approach the flat line E t ~ 0.31 (dot-dash) only after times of many 
thousand lattice units. This means that the mutual equilibration process is associated with 
a very long relaxation time. 

Assuming an exponential relaxation law for the gauge field energy E G (t), 

E G {t) = £ G (oo) +£ expH/T GH ) , (3.10) 

with an adjustable constant Eq, we can determine the relaxation time tgh for the mutual 
equilibration process from a linear fit to \n[(E G (t) — E G (oo))}. 111 Fig- we pl°t hi[(i? G (t) — 
E G (oo)) / Eq(oo)] as a function of time; E G (oo) is taken as E t from Eq. (|3.9| ). The solid line 
corresponds to the actual data while the dotted line represents the linear fit. From its slope 
we extract tgh = 6002 ± 72 lattice units. 

As shown in ( [2.25| ), the rate of energy transfer from the gauge sector to the Higgs sector 
is given by the sum s 2ZE° Ei j° Ei . In a large system with no correlations between the signs 
of E% { and this sum would be a fluctuating quantity with a zero mean value. In our 
calculations, the sum turns out to be fluctuating around a small negative value, which means 
on average there is energy flowing from the gauge sector to the Higgs sector although the 
flow is very slow. Furthermore, from ( [2.5|) and (|2.23|) , we notice that J^ExiJxi i s initially 



proportional to 5 and to R 2 . Both are restricted to small values if we want to work both 
in the weak coupling regime and close to the continuum limit. For larger couplings (i.e. 
for increasing R), keeping 8 small to preserve the proximity to the continuum limit, one 
expects faster mutual equilibration. This is borne out numerically: increasing R by a factor 
10 at fixed S, we found a much shorter relaxation time tgh ~ 100. This is, however, still 
much longer than the time scales for the thermalization of the gauge and the Higgs fields 
individually, which did not change appreciably compared to the case with the smaller value 
for R. 
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IV. DETERMINATION OF THE LYAPUNOV EXPONENTS 



Our main goal in this Section is to determine from the numerical integration of the 
equations of motion the Lyapunov exponents associated with the gauge and the Higgs fields, 
and to relate them to the damping rates obtained in perturbative calculations at finite 



temperature [11 



Let us begin with the discussion of two technical issues. First, when determining the 
Lyapunov exponent from the exponential growth of the distance between two initially nearby 
configurations, one can in principle use different distance measures. However, for all possible 
choices the exponential growth in the distance will be dominated by the eigenmode associated 
with the largest Lyapunov exponent, which is unique to a dynamical system. Unless the 
chosen distance measure accidentally projects out this most unstable mode, different distance 
measures should thus yield the same Lyapunov exponent. (Depending on the power with 
which the classical fields appear in the distance measure the actually measured growth 
rates can be different, but they are always related to the Lyapunov exponent by a simple 
multiplicative constant.) Based on these remarks one expects the two distance measures 
given in ( [2.29|) and (|2.30| ) to yield the same Lyapunov exponent. 



Second, in order to relate results from numerical integration to those from perturbative 
calculations, we need to study the scaling behavior of the Lyapunov exponents with the 
energy (temperature) of the system. For the pure SU(2) gauge system the relation 

X G a « \g 2 E v a (4.1) 
o 



was established ||, with E p = 6Eq/3 denoting the gauge energy per plaquette. In [11 



it was shown by a perturbative calculation with hard thermal loop resummation that, at 
leading order in the coupling constant, in the coupled SU(2) Yang-Mills-Higgs system the 
damping rate of a static gauge boson is the same as that in pure SU(2) gauge theory: 

7G (0) = 0.176 2 T. (4.2) 

For the Higgs damping rate the relation 
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7H (0) = 0.018 g 2 T (4.3) 



was found [|TT|] to hold for small values of the Higgs self coupling A. 

Using the scaling law ( |4.1| ) together with the fact that in weak coupling the energy per 
plaquette E p can be linearly related to the temperature of the system, the authors of p|JT0[ 



found that in the pure SU(2) and SU(3) gauge systems the scaled largest Lyapunov exponent 
coincides with twice the gauge boson damping rate. What should one expect in our case? In 
Ref. |]TTJ it was speculated that here the gauge boson damping rate should be related with 
the Lyapunov exponent extracted from the gauge distance measure ( p.29|) while the Higgs 



boson damping rate could be connected with the Lyapunov exponent extracted from the 
Higgs distance measure (|2.30|) . Since the numerical results did not support this expectation, 
the authors questioned the existence of a simple relationship between damping rates and 
Lyapunov exponents in general, raising the possibility that the coincidence found in pure 
SU(2) and SU(3) gauge theory was a numerical accident. 

Since, however, the largest Lyapunov exponent is related to the most unstable mode in 
the entire coupled system and thus, as discussed above, is picked up by essentially every 
distance measure one can choose, we argue that one should rather expect both measures to 
yield the same Lyapunov exponent. We will show promptly that the numerical computations 
confirm this expectation. If this is the case the Lyapunov exponent can be related with at 



most one of the two damping rates. Since the gluon damping rate (|4.2|) is much larger than 
the Higgs damping rate ( |4.3|) it is clear that if there is such a relation it should be between the 
gluon damping rate and the maximal Lyapunov exponent. To which classically computable 
observable the smaller Higgs damping rate can be related remains at this moment unclear. 
Of course, it could also be that the Lyapunov exponent is related to neither of the damping 
rates. To clarify this issue a careful study of the energy scaling of the Lyapunov exponent 
is indispensible. 

In the upper part of Fig. [| we show that, within the numerical accuracy, the two distance 
measures ( |2.29| ) and ( |2.30| ) indeed yield the same Lyapunov exponent. (Please note the 
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suppressed zero on the vertical axis.) Thus they both measure the largest Lyapunov exponent 
of the entire coupled system. We also see that the Lyapunov exponent shows a slight, 
but characteristic long-time variation which is the same for both distance measures. This 
behaviour will be further discussed below. 

To understand the magnitude and physical meaning of the extracted Lyapunov exponent, 
we studied its scaling behaviour with the energy contained in the system. To this end we 
varied the relative energy contents in the gauge and Higgs subsectors as well as the total 
energy of the system, by changing 8 and R within the permitted ranges. We found that the 
Lyapunov exponent (which reflects the most unstable mode in the entire coupled system) 
does not scale with the total energy of the entire system, but rather with the energy stored in 
the gauge subsector. These scaling tests were performed on the basis of numerical runs over 
a few hundred to a few thousand lattice time units. On these time scales the scaling of the 
largest Lyapunov with the gauge field energy per plaquette was found to be accurate on the 
level of 5 - 10%. (Note that in Fig. || the Lyapunov exponent shows long-time variations over 
a similar range of about 10%.) These findings indicate that the largest Lyapunov exponent 
is associated which the gauge field dynamics, and that the chaoticity in the gauge sector 
dominates the dynamics of the entire system, including the equipartition between potential 
and kinetic energy contributions in the Higgs subsystem. 

Returning to Fig. |4] we can now discuss the magnitude of the Lyapunov exponent scaled 
by the energy per plaquette E p in the gauge field sector. The lower part in Fig. ^ shows 
that for very large times of order ten thousand lattice units the Lyapunov exponent seems 
to converge to a constant value slightly above 0.16. This is close to the value of about 
1/6 which, according to the findings of Ref. |J, is expected (see Eq. (f4.1|)) if the Lyapunov 
exponent is related to the gluon damping rate. The observed long-time variation of the 
Lyapunov exponent with final convergence only after many thousand time steps is apparently 
connected with the already observed long-time variation of the relative energy content of the 
Higgs and gauge field subsectors seen in Fig. ^. As the gauge and Higgs field sectors finally 
reach mutual equilibrium after several ten thousand lattice time units, the scaled Lyapunov 
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exponent Xg/E p settles down to a constant asymptotic value of approximately 1/6. Noting 
that E p = 2T after full equilibration and inserting g — 1, a — 1, we obtain 

4^-0-32, (4-4) 
g 2 Ta 

which together with ( |4.2j ) verifies ( |1.3| ) at the level of 10%. A more accurate verification is 
prohibited by the extremely slow equilibration between the gauge and Higgs subsectors in 
our coupled system. The numerically determined Lyapunov exponent is somewhat on the 
low side; according to Figs. [| and |], the magnitude of this discrepancy is directly related to 
the deviation of the gauge and Higgs field energies from equipartition. 

Let us summarize the findings from this Section: We determined the gauge field Lyapunov 
exponent and verified that for the Yang-Mills-Higgs system it is related to the static gauge 
boson damping rate in the same way as for the pure gauge theory. The quantity Ah, as 
defined in ( |2.30| ) and ( |2.31|) , was shown to give the same Lyapunov exponent and therefore 



to be governed by the gauge degrees of freedom, too. Therefore Ah cannot be related to any 
specific dynamics of the Higgs field, and it should not be called "Higgs Lyapunov exponent" 
as suggested in |TT|. In particular, it is not possible to establish a relation between Ah and 



the Higgs damping rate ( |4.3| ). Unfortunately, we have not been able to come up with an 
alternative observable in the dynamics of classical fields which could be related to 7h- 

V. FURTHER DISCUSSION 

We close this paper with a discussion of our results in comparison with the recent work 



presented in ||1 1|| . For small values of the Higgs self coupling A, Biro and Thoma [11 
found that the Lyapunov exponents Ah and Ac deduced from the Higgs and gauge fields, 
respectively, agree within the numerical accuracy. For very large values of the Higgs self 
coupling, on the other hand, Ah was found to decrease, thereby breaking the coincidence 
between Ah and Ac- For the case of small A we have verified the coincidence of Ah and Ag- It 
has a natural explanation since, as we already discussed, in principle any distance measure 
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between two neighboring trajectories in phase space simply yields the maximal Lyapunov 
exponent of the entire system. 

In view of this argument, the disagreement between Ah and Ag at large A at first seems 
surprising. If, however, the distance measure projects out the eigenmode with the largest 
Lyapunov exponent, it can yield a different value for the growth rate. This is apparently what 
happens at very large values of A for the Higgs distance ( |2.32|) used in fill: As already pointed 
out in [fLl]| , for large A the dynamics of the Higgs amplitude fl2.12j ) becomes completely 
determined by the quartic self- interaction term (|2.19ej ), thereby effectively freezing the Higgs 
field amplitude R. Thus this degree of freedom starts to decouple from the rest of the system, 
and it may happen that the growth of perturbations in the Higgs amplitude is no longer 
influenced by the maximal Lyapunov exponent associated with the gauge field, at least, over 
the time scale used for the determination of the Lyapunov exponent. 

We have not studied here the case of large A, for the following reasons: First, if one 



increases A without decreasing the Higgs amplitude one leaves the perturbative regime ( [2.26Q 



Then the computed Lyapunov exponents can no longer be compared with the damping rates 
from hard thermal loop calculations. Second, for very large A the energy density of the 
quartic interaction term fl2.19e| ) and thus also the total energy density become very large. 
As discussed in Sec. [II B| , in this case one runs into lattice artifacts due to the compactness 



of the magnetic energy on the lattice. 



While we agree with Biro and Thoma [11] about the coincidence of Aq and Ah, we 



disagree with their determination of the temperature. Consequently we also differ in the 
interpretation of their results with respect to their connection with perturbative damping 
rates, believing that their rather negative conclusions are not justified. Biro and Thoma (Tl| 
use the Stefan-Boltzmann law for 10 degrees of freedom 

E tot tt 2 (k B T\ 3 

— = Y kBT [l^) (5 - 1} 

to extract the temperature T from the total energy density E tot /V of the system. Here we 
put in &b, h, and c explicitly (elsewhere we take ks = h = c = 1). As can be seen from 
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the occurrence of Planck's constant in ( |5.1|) , this law makes explicit use of the quantized 
nature of excitation energies in a gas of (massless) particles. As such it can not be applied 
to a classical field theory. The classical Yang-Mills-Higgs fields on the lattice correspond to 
a dynamical system of coupled anharmonic oscillators which in thermal equilibrium and for 
weak coupling obeys the relation 

Ek* 10T 



iV 3 a 3 a 3 



+ o(g,X). (5.2) 



Thus the energy per degree of freedom is directly proportional to the temperature, and not 
to its fourth power. 

Even more crucial is our finding that the whole system equilibrates only after a very 



long time. Since the authors of |TT| did not use the rescaling method for extracting the 
Lyapunov exponents, they ran very soon into the problem of distance saturation due to 
the compactness of the gauge group ||, and thus their analysis was necessarily restricted 
to times which were short compared to the mutual equilibration time between the gauge 
and Higgs sectors. If the Lyapunov exponent is determined on such short time scales, the 
total energy of the system does not tell anything about the temperature of the relevant 
gauge subsector no matter whether it is determined from or ( |5.2| ). For this reason 



Biro and Thoma missed the specific connection of the Lyapunov exponent with the gauge 
subsystem rather than with the entire coupled Yang-Mills-Higgs system. Only an analysis 
like the one presented here can uncover the fact that after very short thermalization times 
one can cleanly define separate temperatures for the gauge and the Higgs subsystems, and 
that the maximal Lyapunov exponent actually scales with the temperature of the gauge 
sector already long before global equilibration between the gauge and Higgs sectors sets in. 

In conclusion, we have studied the high temperature behavior of the coupled Yang- 
Mills-Higgs system in the weak coupling regime close to the continuum limit. We found 
that the coupled system undergoes thermalization in stages: while the gauge and Higgs 
subsectors self-thermalize very rapidly, equipartition of the energy between these two sectors 
takes several orders of magnitude longer. We expect and have qualitative evidence for a 
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connection between the rapid thermalization time scale within each of the two subsectors and 
the magnitude of the largest Lyapunov exponent; to establish this relation quantitatively 
requires, however, a much more detailed and computer-intensive investigation. Since we 
found that the largest Lyapunov exponent is associated with the gauge field and scales 
with the gauge field energy, we conclude that thermalization of both the gauge and Higgs 
subsectors is driven by the chaotic dynamics of classical non-Abelian gauge fields. The Higgs 
subsector thus thermalizes via its coupling to the chaotic gauge fields (gluon heat bath) and 
not due its own intrinsic non-linearity via the selfcoupling A. The scaling of the Lyapunov 
exponent with the gauge field energy was used to verify that, even in the presence of a 
coupling to matter fields, the largest Lyapunov exponent still agrees with twice the damping 
rate of a static gauge boson. 



The analytical results of Thoma and Biro [[n|] from thermal perturbation theory show 
that in the coupled Yang-Mills-Higgs system also the Higgs damping rate is proportional to 
g 2 T and thus, from the general arguments given e.g. in |TIJ, is expected to survive in the 
classical limit. At the moment it remains an open question whether it is possible to identify 
an observable in the dynamical evolution of classical Yang-Mills-Higgs fields which can be 
associated with the Higgs boson damping rate in a similar way as the largest Lyapunov 
exponent is related to the gluon damping rate. 
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APPENDIX A: CALCULATION OF PHASE SPACE FACTORS 



In this Appendix we derive the phase space factors presented in Sec. [HI A| . In general 
it is difficult to calculate the energy distribution over the microscopic degrees of freedom of 
an arbitrary system evolving in time. It is equivalent to solving the equations of motion. 
However, as shown in the main part of the paper, the gauge and Higgs sectors of the system 
equilibrate separately very rapidly (albeit reaching different temperatures). In addition we 
are working in a regime where weak coupling is assumed between the gauge and Higgs sectors 
as well as within each sector. This significantly simplifies the task of calculating the energy 
distributions of different degrees of freedom. 

To determine how the electric energy is distributed over the lattice sites for the thermal- 
ized gauge sector at temperature Tq, we neglect its coupling to the Higgs fields and write 
the partition function for the gauge sector as 

if e i + H n 



II dE x,i M U x,i) exp 



'mag 



- G 



5{D%Ei^5{F{U^) 



(Al) 



where the first 5-function takes into account the constraints imposed by Gauss' law, and 
the second one represents gauge fixing of the form F(U X ^) = 0, which involves link variables 
only. fi(U) is the Haar measure defined on the SU(2) group mainfold. 

For the electric energy on a site, we define its distribution function P e i(e) as 



PM) 



Tr 



Z- 1 J n dE a xi dn(U Xti ) exp 

x,i,a 

x 8{D%E b x J S(F(U^)) S(e-IY1 K,M ,i) 
Z- 1 I ]JdEl i( iLi(U X}i )d\ a x dK5(F(U Xji )) 



x exp 



— + * E Kir':, +i«(e-IE W > ( A2 ) 

x,i,a,b i, a 

where we have introduced Lagrange multipliers A" and k to account for the 5-functions. 
The electric and magnetic energies are defined in ( 2.15|) and ( |2.16| ), respectively. Since the 
electric fields appear only up to quadratic order in (|A2|), the integrations over E xi can be 
performed with the result 
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PM ~ / n dv(U Xti )d\ a x dK8(F(U Xti )) det" 1 / 2 IY-L + i*^. 

Tj rp X^=X 

-nmacr -IG \ - / n ab \a\ 
-7T 2^ \ U x,i A x) 



xo I "x,x' "0,0' 



x exp 



'mag 



l^\ U x a ^i 



x,i,a,b 



2(1 + ikTq) f-i v 



me 



n ^(y^^^y) (1 + iktq) 



-9/2 



x exp 



#mag 



7 <: 2 x S 1 + «T G ^ 
The integration over the Lagrange multipliers A" can be evaluated similarly. After that, the 
whole integral factorizes into integrations over the link variables and the Lagrange multiplier 
k. The integration over the link variables is e-independent and can be factored out as a 
constant. We thus obtain 

+00 



+ ine 



(A3) 



-Pci(e) ~ J dn (1 + ihTq) 3 exp (me) 



—00 

J2 



-3 exp(-e/T G ) , 



(A4) 



which yields the electric phase space factor / e i(e) = e 2 . Three points should be stressed 
here: First, the above result is exact for a pure SU(2) gauge system and approximately valid 
for the weakly coupled SU(2) Higgs system. Second, the apparent temperature agrees with 
the true temperature, i.e. C e i = 1. Finally we mention that e 2 is replaced by e N2 ~ 2 for the 
general case of SU(A r ). 

In deriving the distribution for the magnetic energy on a plaquette, we again neglect the 
coupling to the Higgs fields and write 



Pmag(e) = Z~ 1 Ji[ dIZ,d,,( I ■,..,) exp 



G 



x tar!:, I'-!. ;) 8(F(U X J) 5 (e - 4(1 - ±trU Xo , iojo )) 



(A5) 



Assuming weak coupling within the gauge sector itself, we neglect the nonlinear term in 
the Gauss' law, which reduces D • E = to a linear form: V • E = 0. Then the integration 
over the electric fields decouples from that over the link variables and can be factored out 
as a constant independent of e. Hence we find 
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P mag (e)~ jUMU^) exp (-# mag /T G ) 5{F(U x>i )) S (e - 4(1 - |tr U Xo>iojo )) . (A6) 

x,i 

Now the integration is over the link variables U Xj i, while the energy is expressed in terms of 
the plaquettes U x ^j. In order to carry out the integration, we use a technique described in 
|T8| to transform the integral over link variables into an integral over plaquettes [17 : 



P m ag(e) ~ J II dn(U x ,ij) exp (-H m jT G ) ]J6(U C - 1) 5 (e - 4(1 - \tiU X0M)3o )) , (A7) 



where the 5-function 5(U C — 1) imposes a constraint on each elementary cube c. These 
constraints are the lattice version of the Bianchi identities, which correspond to D • B = 
in the continuum limit. U c is an SU(2) matrix defined in terms of the link variables 
associated with the cube c. The transformation from links to plaquettes corresponds to the 
transformation from vector potentials to magnetic field strengths in the continuum limit. 



For details on the transformation and the exact form of U c , we refer the reader to |18P and 



T7J. Now we consider (|Ay|) in two limiting cases: Tq — > and Tq — ► oo. 

In the low temperature limit, i.e. when the field amplitudes are small, the non-Abelian 
part of the constraints becomes negligible and hence the non-Abelian Bianchi identities 
D • B = reduce to their Abelian equivalents V ■ B = 0. On the lattice, they take the form 
of J2ij£ C Bij = 0, where the sum runs through the six plaquettes forming the surface of a 
cube. These are strong constraints which completely eliminate the longitudinal components 
of the B fields and reduce the effective number of degrees of freedom by 1/3. So in this 
limit, the apparent temperature T a = 2Tq/3. 

In the high temperature limit, we expand the constraints in terms of the irreducible 
representations of the group SU(2) |13| : 

S(U C - 1) = Y,^J + 1) Xj(U c ) , J = 0, |, 1, . . . (A8) 
j 

where Xj{U c ) is the trace of U c in the (2J + l)-dimensional representation. For Tq — > oo, the 
dominant contribution to ( |A7| ) comes from J = 0. Taking J = and neglecting those with 
J 7^ is equivalent to ignoring the constraints imposed by the Bianchi identities and treating 
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the plaquettes as independent variables. Then it is not surprising to find that the apparent 
temperature T a obtained in measuring P mag (e) coincides with the actual temperature Tq. 
While it is easy to do the integrations in the limiting cases, it is difficult to calculate (|A7[) 



for intermediate temperatures. Noting that the major effect of the constraints S(U C — 1) on 
the distribution function is to change its logarithmic slope from the actual temperature Tq 
to the apparent temperature T a , we assume that these constraints can be effectively taken 



into account by replacing Tq by T a in ( A7) and write 



Pmag(e) ~ J II MU x ,ii) exp (-# mag /T a ) 5 (e - 4(1 - ftr U XOiiojo )) . (A9) 

Of course, the validity of this assumption is to be justified and the relation between T a and 
Tq to be established by numerical calculations. 

The evaluation of (|A9|) is straightforward. Making use of the quaternion representation 
U = u° + %T a u a for a plaquette U and the explicit form of the Haar measure on the SU(2) 
group manifold 

dn(U) = -=S (hx{rfU) - l) du° du l du 2 du 3 , (AlO) 
one can easily carry out the integration in (|A9|) and obtain 



-Pmag(c) ~ /mag 

(e)exp (— e/T a ) , (All) 

where / mag (e) = J e(8 — e) . The relation T a = C(Tq) Tq can be determined numerically and 
C(Tq) varies from 2/3 to 1 as Tq increases from zero to infinity ||. 

Note that, for the sake of convenience, we have calculated the distribution of the electric 
energy on a site and that of the magnetic energy on a plaquette. Hence the apparent electric 
and magnetic distributions do not coincide even at low temperature where, in principle, one 
can treat the electric and magnetic energy in exactly the same way and arrive at the same 
distributions. However, as seen in Fig. [I], the different results for T a relate to the same value 
of Tq in the electric and magnetic sector. 



29 



To derive how the Higgs kinetic and potential energies are distributed over the lattice 
sites for the thermalized Higgs sector at temperature Th, we neglect both the gauge-Higgs 
coupling and the Higgs self-coupling. Then it is straightforward to derive 



e 



iWt(e) ~ ^2 exp(-e/T H ) , (A12) 



T' 

i H 



from which we obtain /k in (e) = / po t,2(e) = e- 

In principle, perturbative corrections to the results obtained in (TA4 ), ( All|) , and ( A12|) 



can be calculated. For our case at hand, however, it turns out that the formulae derived 
above are in good agreement with the numerical results. This is not surprising because our 
numerical calculations are performed for weak coupling. 
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FIGURES 

FIG. 1. Short-time equilibration inside the gauge and the Higgs subsectors. Four characteristic 
temperatures are plotted as functions of time: the temperature of the electric gauge sector T c \ 
(solid line), the temperature of the magnetic gauge sector T mag (dashed line), the temperature 
of the kinetic Higgs sector Tk in (dotted line), and the temperature of the quadratic part of the 
potential Higgs sector T poti 2 (dot-dashed line). Parameters: 5 = 0.2, R = 0.2, g = 1, A = 1, and 
a = 1. 

FIG. 2. Long-time equilibration between the gauge and the Higgs subsectors. Five energies are 
plotted as functions of time: the electric energy per gauge field degree of freedom E e \ (solid line) , 
the magnetic energy per gauge field degree of freedom E ma _ g (dotted line), the total gauge field 
energy Eq = E e \ + E mag per gauge field degree of freedom (short dash), the Higgs field energy 
per Higgs field degree of freedom E-r (long dash), and the total energy per degree of freedom Et 
(dot-dashed line). Parameters as in Fig. |l[ 

FIG. 3. Determination of the relaxation time tgh associated with the long-time equilibration 
between the gauge and the Higgs subsectors. The solid represents the numerical data while the 
dotted line is a linear fit to the data. 

FIG. 4. Upper graph: coincidence of Aq and Ah- Lower graph: convergence of the gauge field 
Lyapunov exponent. E p is the gauge field energy per plaquette. The fields were rescaled after each 
period At = 30. 
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